Fast and Elegant Numerical Linear Algebra Using the RcppEigen Package
نویسندگان
چکیده
The RcppEigen package provides access from R (R Core Team 2015a) to the Eigen (Guennebaud, Jacob et al. 2012) C++ template library for numerical linear algebra. Rcpp (Eddelbuettel et al. 2017; Eddelbuettel 2013) classes and specializations of the C++ templated functions as and wrap from Rcpp provide the “glue” for passing objects from R to C++ and back. Several introductory examples are presented. This is followed by an in-depth discussion of various available approaches for solving least-squares problems, including rank-revealing methods, concluding with an empirical run-time comparison. Last but not least, sparse matrix methods are discussed.
منابع مشابه
A Rational Approach to Portable High Performance: The Basic Linear Algebra Instruction Set (BLAIS) and the Fixed Algorithm Size Template (FAST) Library
We introduce a collection of high performance kernels for basic linear algebra. The kernels encapsulate small xed size computations in order to provide building blocks for numerical libraries in C++. The sizes are templated parameters of the kernels, so they can be easily congured to a speci c architecture for portability. In this way the BLAIS delivers the power of such code generation systems...
متن کاملComputation of a canonical form for linear differential-algebraic equations
This paper describes how a commonly used canonical form for linear differential-algebraic equations can be computed using numerical software from the linear algebra package LAPACK. This makes it possible to automate for example observer construction and parameter estimation in linear models generated by a modeling language like Modelica.
متن کاملA Rational Approach to Portable High Performance: The Basic Linear Algebra
We introduce a collection of high performance kernels for basic linear algebra. The kernels encapsulate small xed size computations in order to provide building blocks for numerical libraries in C++. The sizes are templated parameters of the kernels, so they can be easily congured to a speci c architecture for portability. In this way the BLAIS delivers the power of such code generation systems...
متن کاملA Rational Approach to Portable HighPerformance : The Basic Linear
We introduce a collection of high performance kernels for basic linear algebra. The kernels encapsulate small xed size computations in order to provide building blocks for numerical libraries in C++. The sizes are templated parameters of the kernels, so they can be easily con-gured to a speciic architecture for portability. In this way the BLAIS delivers the power of such code generation system...
متن کاملClassical Wavelet Transforms over Finite Fields
This article introduces a systematic study for computational aspects of classical wavelet transforms over finite fields using tools from computational harmonic analysis and also theoretical linear algebra. We present a concrete formulation for the Frobenius norm of the classical wavelet transforms over finite fields. It is shown that each vector defined over a finite field can be represented as...
متن کامل